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A directly measurable correlation length may be defined for systems having a two-step relaxation, 
based on the geometric properties of density profile that remains after averaging out the fast motion. 
We argue that the length diverges if and when the slow timescale diverges, whatever the microscopic 
mechanism at the origin of the slowing down. Measuring the length amounts to determining ex- 
plicitly the complexity from the observed particle configurations. One may compute in the same 
way the Renyi complexities K q , their relative behavior for different q characterizes the mechanism 
underlying the transition. In particular, the 'Random First Order' scenario predicts that in the glass 
phase K q = for q > x, and K q > for q < x, with x the Parisi parameter. The hypothesis of a 
nonequilibrium effective temperature may also be directly tested directly from configurations. 

PACS numbers: 64.70.Q-; 61.43. Fs; 63.50.Lm; 64.70.P- 

I. INTRODUCTION 

A solid is an arrangement of particles with permanent density modulations. While commonplace, it is nevertheless 
remarkable that a system of soft particles may condense into a crystalline state, since due to thermal activation, 
no particle is strictly confined to a specific region of space. This notwithstanding, permanent density modulations 
persist forever in an infinite sample. The reason for this is that in order for spontaneous thermal fluctuations to 
completely erase the density modulations in the solid phase, rearrangements involving a divergent number of particles 
are required, and this in turn requires infinite free-energies, and as such is an infinitely rare event. 

Glasses also appear solid, in the sense that they have long-lived density modulations, although it is still not known 
whether these are truly permanent as for crystals. Similarly, relaxation in some supercooled liquids occurs on a 
timescale that appears to diverge faster than Arrhenius as the temperature is lowered. Both these features of slow 
relaxation, if they hold rigorously even for a system of soft particles, may only be explained if there is an ever increasing 
number of particles having correlated evolutions, implying that there is a growing coherence length [I] as a function 
of decreasing temperature. The central puzzle of glasses is that there does not appear to be any recognizable form 
of spatial order - on the contrary, the configurations seem definitely liquid-like, so it is not clear from whence they 
derive their solidity. 

There are several possible resolutions to this: the simplest is to conclude that glasses are not truly solids. That is, 
their relaxation would become purely Arrhenius at sufficiently low temperatures, as is expected for a system comprised 
of independent parts. In this case, an amorphous solid would strictly speaking exist only in the trivial case of hard 
particles immobilized by contact with others, for which all activation is absent. Another possibility is that (at least 
some) glasses indeed have a truly super- Arrhenius behavior as T — > or even as T — > Tk, for some finite temperature 
Tk, termed the "Kauzmann temperature" [Bl. This assumption is implicit in several theories of fragile glasses. For 
example, density functional theory expresses the liquid properties in terms of a free energy functional of the the 
density, and exhibits, for low temperatures and high pressures, not only periodic 'crystalline' configurations, but also 
non-periodic, amorphous states [5]. A related, more complete version of this is the 'Random First Order' scenario 0[7j, 
which posits an 'ideal glass' state to which the system would equilibrate given an infinitely slow annealing from high 
temperature to T < Tr. Such an ideal state would presumably be characterized by permanent non-periodic spatial 
density-modulations. In order for the solids postulated in these theories to exist in finite dimensions beyond the 
mean-field approximation, they require coherence at diverging lengths [2]. This leads to a central question: Can this 
cooperativity be observed directly from configurations? 

In this paper we shall discuss a family of lengths [5] that, we argue, should diverge in any system having a timescale 
that diverges at some temperature (or equivalently, shows a super- Arrhenius behavior as T — > 0). The lengths are 
geometric in nature, in the sense that they do not require information on either microscopic dynamics, the history of 
the sample, or even whether the system is in or out of equilibrium. The various scenarios for the glass transition are 
distinguished by the relative manner in which the different lengths behave. 

The reasoning proceeds as follows. First, we recognize that order concerns the time-averaged density, which is 
insensitive to fast excitations such as vibrations, particle-vacancy excitations, or isolated spin flips. In fact, this is 
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true as well for crystals, which are only truly periodic once such excitations are averaged out [5]. In a non-crystalline 
system, these average density profiles (which we shall term 'states') are not periodic, so a multitude of possibilities 
for the density profiles are obtained simply by translations. This leads naturally to the question: How many possible 
states are there? The logarithm of this number is the complexity (we reserve the word 'entropy' for quantities based 
on counting individual configurations, and 'complexity' for those associated to classes of configurations). We note 
that care must be taken here, because density distributions are continuous objects, so that counting them is much 
like counting the number of trajectories of a dynamical system. We shall exploit this analogy. 

At this point we make the following observation: if we wish to concentrate on a truly solid phase, in which a 
density profile lasts forever, the number of different possible density profiles cannot be exponential in the volume, or, 
to be more precise, the complexity cannot be extensive. This is because if there are exponentially many metastable 
states in competition, a simple nucleation argument suggests that they cannot all have infinite lifetime; starting from 
any particular state, some other state will nucleate in finite time. (This is of course only true for short range, finite 
interactions in finite dimensions.) As we shall argue, the reason is that if the number of possible finite-size motifs the 
system may nucleate is exponential in their volume, the increase in probability of nucleating a new motif overcomes 
the cost, in terms of probability, of creating an interface. 

We now consider a sample of very large size V, and 'patches' of volume V <C V belonging to it. The fact that the 
number of density profiles in a solid is subexponential immediately implies that patches of configurations of volume 
V repeat more often than exponentially rarely in V, so that their number Af(V) obeys ln(Af(V)) < 0(V). On 
the contrary, within the liquid phase, or in a defective (nonequilibrium) glass, the system breaks into uncorrclatcd 
regions, and we have, for large V, \n(JV(V)) ~ VK\, with K\ being the complexity per unit volume in the large V 
limit (limy-yoo KY), with the statistics taken over patches belonging to the (much larger) volume V. 

The fact that patterns of all sizes repeat less often than exponentially, so that there is no (extensive) complexity, 
directly implies the existence of a divergent correlation length. We may see this most easily as follows. Suppose we 
have a region of volume V with a configuration A. To what extent does A determine the configuration (say, B) of a 
neighboring region, also of size VI The typical number of repeats of configurations of volume V (or rather, its typical 
logarithm) is \nAf(A) ~ \nJ\f(B) ~ — K\ . Clearly, if A and B appeared independently, the repetition probability 
of the pair A + B would be ]nAf(A + B) = \nU{A) + ha.N'(B) ~ -2VKY, while in general it will be some number 
lnAf(A + B) ~ —VKi V . The patches will then not be independent if Kf v < 2KY , or, in other words, if the entropy 
of patches of size V is subextensive. In that case, the knowledge of the configuration in a volume V gives information 
on that of an adjoining patch. Note that we have assumed here that for large V, the most likely patch entropy (based 
on the logarithm of the frequency of appearance of a patch) is, to leading order in V, independent of the shape. If we 
have an extensive patch entropy K\ , and we assume that for large volumes the system may be considered as composed 
of a number of independent components of volume V\ ~ , then Af ~ e VKl ~ e v / Vl , and we identify a correlation 
volume as V\ ~ 1/K\. 

A system with subextensive patch entropy K\ will then have an infinite length - in this sense, it is ordered. Now, 
let us start from such an ordered system, and fracture it into pieces of typical size V c . Let us now reassemble these 
pieces at random to form a new system of large size. In this case, we expect patches smaller than V c will be found 
often, and those larger than V c will be exponentially rare, as in Fig. [T] This leads us then to the definition of the 
following types of length for an amorphous system: 

i_ 

• Complexity length: l\ — K 1 D , measures the frequency with which a typical patch repeats. 

• Renyi lengths: In order to distinguish different glass transition scenarios, it is useful to consider patches 
that repeat more often or less often than is typical. This leads to the definition of lengths based on the Renyi 
complexities, or, equivalently, the large deviation function for patch repetition frequencies. 

• Crossover length: patches with volumes smaller than V -C V c repeat often, while larger ones V c <C V are 
exponentially rare. The crossover length is i ~ V C D ■ 

A related measure is the following: suppose the complexity of a patch of volume V is given, with finite V 
corrections, by: VK\ — A(V)k[ 1 " > + where A(V) is subleading in V, for example A(V) ~ V. The next to 
leading term A(V)k[ 1 " > describes correlations, and is hence referred to as the 'synergy' [8]. K± might diverge 
at the critical point jMj. If A(V) is not O(l) in V, provides a new length. 

These measures may be understood easily in the case of a microcrystalline system, as sketched in Figure [TJ However, 
one of the central points of our proposal is that it allows us to define a correlation length for systems whose ground 
states are not only non-periodic, but even conventionally thought of as amorphous, in the sense that their scattering 
patterns have no Bragg peaks. In what follows, we shall discuss in detail how the counting and identification of patters 
may be made, as well as the implications to current glass scenarios. 



FIG. 1: Two patterns, one smaller and one larger than the crystallite length. The small pattern will recur with high frequency, 
but the large one will be rare, because it requires recurrence of the crystallite boundaries within a patch, which is random 
event. Patches that are very large with respect to the crystallite volume occur with multiplicity exponential in the number of 
defects per volume, of the order of the number of microcrystallites contained. 




FIG. 2: Looking for a patch that is the same, up to a certain precision, and allowing for rotations. 



II. SCRAMBLED SEQUENCES AND IMPERFECT TILINGS 



A. Perfect deterministic sequences 



It is important to note that sub-extensive complexity of patches does not imply that the infinite system is either 
periodic or quasiperiodic, as characterized by a diffraction spectrum consisting of S function Bragg peaks. To make 
this clear, we review a particularly well understood class of systems, namely those which possess 'recursive symmetry'. 
For ease, we shall consider one-dimensional sequences of two symbols (1 and 0) which may be generated by means of 
a substitution rule, wherein a symbol is replaced by a specific word (composed of 0's and l's) everywhere it occurs, 
and then iterating this process. Generalization to higher dimensions, for which building blocks are subdivided, is 
immediate. We consider three prototypes which characterize three types of spatial structure: 

Rule 1 : 1 ->■ 10 ->■ 10 (1) 
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Starting from (say) 0, we generate sequences of length 2™ by iterating this rule n times; the result is clearly periodic: 



0101010101010101010101010101 



The second possibility is illustrated by 

Rule 2 : 1 -> 10 -> 1 (2) 

If we start from as before, we generate strings of length F n , where F n is the n th Fibonacci number: 

1011010110110101101011011010110110 

This rule produces the Fibonacci sequence, which is quasiperiodic, with a discrete Fourier transform consisting of a 
dense set of S -functions. Finally, we consider Rule 3: 

11 -> 1110 

10 1101 

01 ->■ 0010 

00 0001 

Iterating this rule n times produces a string of length 2™ : 



1110110111100010000111011110110111101101111000101110110100011101 (3) 



This sequence, known as the Rudin-Shapiro sequence, has a discrete Fourier transform with flat modulus [10h13| : and 
Rule 3': 

1 -> 10 -> 01 (4) 

which produces: 



01101001100101101001011001101001 



(5) 



This sequence, known as the Thue-Morse sequence, also has no ^-functions in its Fourier transform |lll I12j The 
Rudin-Shapiro and the Thue-Morse sequences are examples of "Non Pisot sequences"[10 . 
Thus, the three examples have different kinds of order: 

1. Periodic: Fourier spectrum peaked at the period and its harmonics. Degeneracy = 2, (entropy = In 2) (and, 
more generally the logarithm of the cell size). 

2. Quasiperiodic: Fourier spectrum is dense set of <5 -functions. Two sequences taken at random have an overlap 
Q distributed according to P(Q) consisting of a 5 function and a uniform tail, as in Fig. [3] the Parisi function 
P(Q) is non-trivial. The number of subsequences of length £ (that is, with I symbols) appearing in a large 
sequence is proportional to £, (actually I + 1) as is easy to see from the rules above. All patches of size £ repeat 
within a distance of 0(£). The probability that any two subsequences taken at random have a Hamming distance 
per unit length < e is proportional to e, and independent of their size £. 

3. Non-Pisot sequence: For the example of the Rudin-Shairo sequence, the modulus of the discrete Fourier 
transform is flat, independent of u. The Parisi function P(Q) is peaked: two (asymptotically large) subsequences 
taken at random from the infinite sequence have an overlap that is, with high probability the smallest value 
possible (one-half, in this case). Here again, the number of subsequences of length £ is also oc £, and all patches 
of size £ repeat within a distance of 0{£). 

In going from 1-3, some of the conventional (scattering) manifestations of order are lost [? ], despite the fact that 
the entropies of sequences of types 2 and 3 are still logarithmic, and patches of size £ repeat much more often than 2 , as 
they typically would for a random system. To understand this frequency of repetition, note that a subsequence of size 
£ is completely included within the 'descendants' of a single site after k oc log(() substitution steps; the subsequence 
repeats in the descendant of a neighboring ancestor - its 'cousin /c-times removed'. This argument is applicable to all 
systems obtained by substitutions jTHJ, [T7] . 

We conclude that both the crossover length £ c as well as the complexity length l\ diverge in all these deterministic 
sequences - even though there may be no Bragg peaks in the diffraction spectrum, so that by the conventional measure 
of scattering such systems would be considered amorphous. 
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average overlap 



FIG. 3: Probability of overlap of two Fibonacci sequences P(q 



(a) 



(b) 



FIG. 4: (a) A set of Wang tiles with matching rules indicated by numbered edges. The rule for laying the tiles is that only 
edges of the same number may abut, (b) States obtained from a Monte-Carlo simulation of a Wang tiling quenched from high 
to low temperature are left with only 'hard' defects (i.e. mismatches of adjacent sides); as indicated by the red lines. 

These considerations may be generalized to higher dimensions as well. Perhaps the best known 2D quasiperiodic 
tiling is the Penrose tiling, but a system with similar properties is a Wang tiling |18j . In constructing a Wang tiling 
one is required to follow matching rules. One set of tiles with their associated matching rules is shown in Figure ^sl): 
every interface between tiles has to have the same number on both tiles, a two-dimensional domino. 

An example of a system that is statistically isotropic is a tiling of identical triangles 19! called the 'Pinwheel' tiling 
(Fig [6]) . The Pinwheel tiling can be constructed by recursive substitution rules, so patterns of size I repeat (modulo 
rotation) every V ~ I in a perfect tiling, as for the ID sequences discussed above. 

Significantly, it is known that any packing which can be generated by a recursive substitution rule admits a set of 
matching rules which only allows tilings consistent with the recursive symmetry [20] ; the set of matching rules for the 
Pinwheel tiling is known explicitly [19] . 

Two-dimensional tilings of the Non-Pisot type have also been constructed [21] . This means that an experimental 
system having this kind of order would go unnoticed, if the only way to detect quasicrystals order is a diffraction 
pattern! The strategy we shall outline in the following sections may be used to detect this form of order. 

B. Imperfect sequences: Finite lengths 

The examples of the previous subsection are all systems with infinite correlation lengths. What does an imperfect 
quasiperiodic system, or an imperfect Non-Pisot system look like? The question begins to be tricky, as there is 
no notion of periodicity in the latter case. In analogy to a system which has not achieved the ideal glass state, 
let us consider a sequence constructed by concatenating unrelated subsequences of varying lengths n, drawn from a 
distribution P(n), with average n. For example, the scrambled Thue-Morse sequence shown below: 



01101 0110100101 0011001 1001101001 



(6) 



is the analogue of a polycrystallite, but with micro/Non-Pisot pieces. It is easy to see that strings of length I -C n 
will repeat often, because even in independent subsequences the same string will appear. However, if t 3> n, the 



G 



<v(t)> 



FIG. 5: The average distance T>(£) between words of length I in a perfect Thue-Morse sequence of length 1000 (blue) and a 
sequence of length 1000 composed by concatenating strings of average length 50 randomly drawn from a perfect Thue-Morse 
sequence (red). 

likelihood of overlap is exponentially small, by the central limit theorem. This is seen in Figure [5] where the average 
distance between strings of length I is plotted as a function of £ for a perfect and a randomly concatenated Thue-Morse 
sequence. The two lengths l c and £\ are then estimates of the coherence length n. 





FIG. 6: The Pinwheel tiling, with infinite correlation length. 



Consider next the example of the Wang tilings, which has the virtue that it can be expressed as a lattice model of 
spins [22 (taking 16 values in the example of [12]), assigning an energy corresponding to a mismatched edge - that is 
to say, the matching rules may be enforced by a suitable Hamiltonian. Consequently, such a system can be simulated 
at finite temperature (See Figure |4jb)) to measure the growth of the coherence length. At finite temperature, the 
Wang tile system has energetic point defects - local violations of the matching rules. Starting from a high temperature 
state and annealing down to low temperatures, many easy-to-cure 'benign' defects disappear, leaving the system in a 
state with isolated 'hard' defects (Figure |4jb)) that can only be removed by large-scale rearrangements of the order 
of the interdefect distance. Such behavior is also present in the Penrose tilings, where there are benign matching 
rule violations due to local tile flips, and non-trivial configurations of the phason field |23] which require large-scale 
rearrangements to eradicate. In this case, as we have verified numerically for the Wang tile system, patches repeat 
often when much smaller than the distance between serious defects, and rarely when they are larger. 
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III. PARTICLE SYSTEMS: DEFINING AND COUNTING DENSITY PROFILES 

A. Density profiles 

Particle systems are more subtle to compare than number sequences, since they have continuous positions and 
thermal vibrations. We will now present a method that allows us to distinguish and count configurations in principle. 
Recalling our motivation from glasses, we proceed in two steps: first we average out rapid fluctuations to obtain a 
continuous density profile (a state), and then we count the number of profiles. 

Consider first a snapshot of a crystal at finite temperature. It is clear that, strictly speaking, no two patches of 
the crystal coincide exactly, because of thermal fluctuations in the atomic positions. Therefore, we consider a density 
profile averaged over a long (in principle infinite) time: 



1 f Tp 

p(r) = lim — / dt p(r, t) 



(7) 



with 



p(r,f)=X;*(r-r B (t)) 



(8) 



density autocorrelation (t-t ) 



t-t 



POO 



FIG. 7: A sketch of the situation in a crystal. At finite temperatures, particles are not in periodic positions. However, the 
average density, defined by performing an infinite time average, is indeed periodic. In this case, we also have the possibility 
of averaging over spatial repetitions, a procedure automatically performed by the Fourier transform. The density-density 
correlation function has a plateau in time, which lasts forever if there is positional order. 



In the long-time limit (in this paper, this always means taken after the thermodynamic limit) , the density tends to 
a smooth function with a well-defined crystal modulation, as sketched in Figure [7] The fact that there are permanent 
density modulations is reflected in the fact that autocorrelation function 



c(t,t') = v 



dr (p(r,t)-p) (p(r,t')-p) 



(9) 



has a plateau qEA as t — t' — > oo: 



i J dr (p(r) - pf = i / ,11, II' C U.I') 



qEA > 



(10) 



where 



p = ~ J dr p{r) 



and ~ denotes in the large r p limit. This defines the Edwards- Anderson parameter qEA- 

An ideal amorphous solid should have permanent, non-periodic density modulations. By this we mean that, if 
we define a density profile as in equation ([7]), we obtain a function that may be neither constant nor periodic nor 
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quasiperiodic. Such solids might or might not strictly speaking exist. What we know in practice is that a supercooled 
liquid has a relaxation time r Q , much larger than the microscopic time r m , during which an amorphous density profile 
remains stable, so that the autocorrelation function ([9]) has a plateau as in Figure [8] Only if r a — > oo clo we have a 
true amorphous solid. 



C(t-t') 



V 




t-t' 



FIG. 8: A sketch of the autocorrelation function with a two-step relaxation, plotted in logarithmic scale for the times. There 
is a rapid decrease into the a plateau, whose value is by definition the Edwards- Anderson parameter, followed by a fall to zero 
at a longer time scale r a . The slow timescale r a for an equilibrium system depends on the temperature, and grows as the 
glassy phase is approached. Out of equilibrium, r a grows with the 'waiting' time after preparation, and saturates if and when 
equilibrium is reached. 



For finite r Q , we may define a density profile for a supercooled liquid configuration as in equation Q, but with the 
average taken over a finite time r p such that r m <C t p <C t q . This brings about an essential limitation on the precision 
with which both p(r) and r are defined. The density resolution is of the order of, but inevitably not better than: 



<Tp) 



dvdt(p 2 (v,t)-p 2 (v)) 



2 






[hi 



dt dt' C{t,t') - q EA 



(11) 



The precise manner in which the resolution e(r p ) goes to zero depends on the way that [C(i — t') — qEA] goes to zero 
on approaching the plateau: if it does so in an integrable form, then e(r p ) ~ (r p )~i. 

The resolution of r is in turn limited by the fact that the density is generated by time-dependent 5 functions , 
corresponding to particle positions. Suppose we divide the total volume V in cubes of side Ar. A particle travelling 
with velocity v will visit vr p /Ar cubes in time r p . Estimating v ~ ^/fc^T, and considering that there are N particles, 
a total of order N\/ksTT p / Ar are visited. If we want this to be larger than the total number V/[Ar] D of cubes, we 
get the condition: 



Ar(r p ) ~ ^fk^TTpp 



(12) 



For large r p , Ar(r p ) is exceedingly small, and we shall assume throughout that we have discretized space or smoothed 
the average density functions up to a distance of this order. In what follows we shall, in addition, sample p at even 
larger spatial intervals, but this we shall do after having discretized space at the scale (12 1. 



In realistic situations, the above time-averaging procedure may be impractical. Two alternatives to this procedure 
are 

• Identifying positions modulo some fixed tolerance that we know a priori to be appropriate, for example a fraction 
of the inter-particle spacing. 

• Instead of comparing the averaged configurations, we may omit the time-averaging and compare the inherent 
structures of configurations, that is, the positions reached after a quench to T = 0. Using inherent structures, 
however, is not free of artifacts and ambiguities: two inherent structures that differ in buckling modes that flip-flop 
on a timescale < r Q should be treated as being the same, while a naive counting will treat them as different. If the 
density of such localized flip-flops is finite, one gets an entropy that may be absent in a correctly time-averaged 
procedure. 
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B. Correlation lengths 



The preceding procedure yields a density profile that has averaged out - to the extent possible - all rapid degrees 
of freedom, such as particle- vacancy pairs, phonons, particle rattlings. Our problem now is how to count the resulting 
continuous profiles, or, alternatively, how to identify congruent patches. This situation is already familiar in a different 
context: it arises when we wish to determine the entropy of trajectories of a dynamical system, in which setting it has 
been well discussed both in principle and in practice. In the case of dynamical systems, we first discretize phase-space 
in cubes of size e and time in intervals r. We consider pieces of trajectories d time-intervals long, and define the 
Kolmogorov entropy on the basis of the probability P e (*ij *d) that at time r the trajectory lies in the cube i%, at 
time 2t it lies in the cube %2, • and at time dr it lies in cube id, as |24) : 

Kx = - lira Una lim — V" P t (i 1: ...,i d )\nP e (i 1 , ...,i d ) (13) 

T->0c->0ii->oo TO z — ' 

More generally, we may define all the Renyi complexities: 

K„ = - lim lim lim — r In [ V P e (h, i d ) q I (14) 

Note that d goes to infinity before e and r go to zero, a fact that shall turn out to be important. 

In practice, one may proceed as proposed by Grassberger and Procaccia |26j : we consider a very long sequence of 
M times, and for every patch of length d labeled by its starting point i (Fig. [To] ) , we count the number of patches 

( e ) th & t coincide with it - up to precision e. The entropies are estimated by 

K a ~ - lim lim lim — ^ In J — VVW 9-15 \ (15) 



and in particular: 



K l ~- lim lim lim — <^ — > lnfn ?>] ) (16) 



This way of putting the problem is practical, because we may extrapolate the results from finite e (for which we have 



many coincidences) to smaller e. This is illustrated in Fig. 12 (a sketch) 



The fact that an entropy may be defined and computed may be understood if we assume that [24] the dependence 
on e (for small e) is subdominant in d: 

y^[n^] 9 - e Td ^+i e ^(9) = e rdK q+1 +4>(q) hie ^ 
i 

for some 4>(q). As we shall discuss later, this means that d goes to infinity keeping | lne|/e? small. 

We may easily adapt this procedure to count density profiles. We consider average profiles p(r) in D-dimensional 



space (Fig. 11 ). We discretize space in some manner, with a typical grid spacing A playing the role of the time-interval 
t in a dynamic system [25 , and we also discretize p in intervals of size e. Here again, we may proceed as follows. 
We consider a large sample containing M different patches ai, ...,aM each of volume V. We look for the number of 
coincidences between a patch and all others, centered about any site and allowing for rotations. Two patches coincide 
if the sum of all the differences in the value of p between corresponding sites is smaller than e. Denoting by (e) 
the number of patches that coincide in this way with the patch a, we have: 



if„~- lim lim lim -, ^— In J — V (e)} ( - q - 1) 1 

i A^oe^ov^oo (q- l)V 1M^ 1 " WJ J 



(18) 



and in particular: 



K t ~ - lim lim lim - \ — Vlnk( y) (e)l > (19) 
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FIG. 9: Counting trajectories in a discretized space and time. 




FIG. 10: Discretization of trajectories. A sequence is repeated three times at this level of precision, but each one will be 
counted as distinct at a smaller e. 



Here and in what follows we shall for (conceptual, though not practical) simplicity keep the limits A — > and e — > 0, 
although they are probably unnecessary provided that one takes V — > oo. 

The Renyi complexities are at most extensive (K q = 0(1) as V — > oo), or subextensive, and then K q — > 0. Writing 
equation ( 18 ) in the form 



ivk 9+1 „ J_y [n (v) {e)]q 

M ^ l 



M ^ 



(20) 





FIG. 11: The analogue of a trajectory is the density profile. The density plays the role of the position in a dynamical system, 
and space the role of time. 
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(where the symbol '~' means ' in the limit A — >• f-)0 V —> oo') suggests that a large deviation function be 
defined for the variable BY = (the sign and the factor B entering the definition of Y are purely conventional) . 

Introducing the large deviation function G through: 



V\BY = - 



ln ii a 
V 



,v G(Y) 



La 

dq e ^W-K^\ 



(21) 

Thus, the Renyi complexities may be thought of as (a form of) a Legendre transform of the large-deviation function 
G(Y). We may evaluate the integral in (21 1 by saddle-point 



G(Y) = q*[f3Y - K g , 



+U 



with BY 



dq 



[qK, 



9+1. 



(22) 



Before proceeding, let us be more precise about the dependence on e. If we assume, as in (17), that (20) may be 
written in the limit of small but non zero (A, e) and large but not infinite V as: 



(23) 



for small e. Note that there might be a correction intermediate in order V between the leading one and the one 
dependent on e. If we work under the condition 



ln< 



VK, 



< 1 



e » e 



-VK q+1 /</> 



9+1 



we may refine equation (22) to read: 

\a.V{Y) = VA G{Y) +lne G D (Y) with G Q = (f>(q*(Y)) 



(24) 



(25) 



with q*(Y) given in (22). From this equation we read the approach, as e e Kq+lV >> 0, of the large-deviation function 
G to its large-]/ limit. Let us stress that we do not know what exactly the dependence of entropy on e might be like 
in a true glassy system, and this is one of the interesting questions an experimental/numerical investigation should 
answer. 



C. Direct measurement of the complexity 

The (patch, or block) complexity Yj V , is defined as the log of the number of different kinds of states for a given patch 
size V, per unit volume. (Here, and in what follows, £ without supraindex stands for the large V limit). Because 
each state a is repeated n a times, we have to divide by n a to count each type once 

E(Y)= lim lim lim -m/ ^"^ ! = BY + G(Y) and BY = (26) 
A^oe^ov^cvV [ n(e) J V 



From equations (26) and (21) we obtain the useful expression: 

e -(,-i)VJf a = f d Y e s ( y )"^ y (27) 



In practice, one would work at different values of A and V , and make histograms of the probability that patterns 
repeat with frequency n(e) at given precision e. Plotting in terms of the normalized logarithm of this frequency 

BY = In (^-^f^J one should obtain a set of curves that scale correctly as E(e, V) — ^(Y) + ^ G (Y). 

The content of the functions G(Y), S(F) and K q 's is the same. They encode important information to distinguish 
different scenarios for the amorphous state. 
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FIG. 12: A sketch of the plots of the logarithm of patch frequency versus precision e, for different lengths d (see Ref [26] )• The 
Grassberger-Procaccia procedure takes advantage of the scaling with e to extrapolate results from relatively large e, to small e, 
a limit in which no coincidences would in practice be found. 

IV. LONG RANGE ORDER 
A. Amorphous order vs. Weak-periodicity 

The property of weak periodicity, introduced by Aubry, is denned as follows [28]: for patches of size V, there is 
a distance R(V, e) along the sample such that in an infinite sample every patch will recur - up to precision e - at a 
distance closer than R(V, e). Note that nothing is said about the dependence of R{V, e) with V. What this definition 
excludes is the possibility that for any given type of patch a, there might be, in an infinite sample, arbitrarily large 
a-free regions, something that would clearly happen in a random sample. Aubry [28] showed that, under reasonable 
assumptions, the ground-state configuration of a particle system with short range interactions is weakly periodic. 

In this paper, we are defining a length l\ for amorphous order also based on patch repetition which is, although 
conceptually related, not equivalent. Central to our argument is the necessity of time- averaging, which allows us 
to compare average-density profiles. We are measuring complexities - based on the multiplicities of these averaged 
profiles - rather than entropies, which are related to the multiplicities of configurations. Time averaging allows us to 
work at non-zero temperature, a situation in which weak periodicity of configurations necessarily breaks down. 

The definition we follow in this paper of perfect order, implying an infinite correlation length, is a situation in which 
almost every patch recurs often, the recurrence distance being less than exponential in the patch's size. This is related 
to the vanishing of the complexity, and it does not exclude the possibility that there will be arbitrarily large regions 
where some patch is absent. 

Is it possible that we could use the notion of weak periodicity applied, rather than to configurations, to time- 
averaged density profiles? The answer is that, even in that case, the two conditions (vanishing complexity and weak 
periodicity) are distinct. To understand this, consider first the case of a quasicrystal at low but non-zero temperature, 
for example the Wang-tile system discussed above. In equilibrium, defects will appear and disappear continuously, just 
as in a crystal. These 'benign' defects will be averaged out in our procedure, so we may safely ignore them. However, 
there are also 'serious' defects having low energies that cannot be eliminated with a finite number of rearrangements 
(the best known being the 'decapods' of the Penrose tilings). These are long-lived defects, and will therefore show up 
in a time-averaged density function. Moreover, they have, in equilibrium, finite probability if their energy is finite; as 
is the case in the Wang-tile and Penrose system above. 

Denote by f a the free energy density of the metastable state a (measured with respect to the one f Q of a perfectly 
ordered system) defined as having a fixed set of 'serious' defects, and including all rapid thermal fluctuations. f a 
is computed from the partial Gibbs-Boltzmann distribution, keeping 'serious' defects fixed. Denote Af = e VE '^ the 
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number of different such metastable states of free energy density / for a sample of total volume V. We may write the 
partition function at low temperature as: 

Z = e-W. [dfe v ™-M (28) 



The integral may be evaluated for large V. If > 4 a saddle point with / > dominates, there is a finite 

aj /=o 

density of 'serious' defects in equilibrium. (As we shall see below, a finite density of defects brings the system back 
to the liquid phase, in this case because defects are so close that they are may be mutually annihilated with finite 

rearrangements, and have no longer infinite life). If, on the other hand ^f 

order is not destroyed because the density of defects tends to zero, although the total number present may still be 
large. In the latter case we have an ordered system with infinite complexity-related correlation length, but no weak 
periodicity (even for a time-averaged configuration) since the 'serious' defects spoil it. 

Although we have chosen to discuss this situation in terms of the Wang tile system, it is also relevant for the (more 
hypothetical) 'ideal glass' scenario at finite temperature, as we shall see. 



< A, the system is frozen in / ~ 0, 

/=o 



B. Infinite lifetime and diverging lengths. 



We shall now argue that if a system of particles with short-range smooth interactions has a state with permanent 
density modulations, then these density profiles necessarily have subextensive complexity, and hence there is, by our 
definition, long-range amorphous order. This is also true for the zero-temperature limit of a system with super- 
Arrhenius timescale behavior. Our reasoning is not rigorous, but we believe it may be made so. The argument is 
based on an elementary observation: it is not possible to have exponentially many metastable states all of which have 
a divergent lifetime. 

A situation that arises in supercooled liquids is that a system has many options of phases for nucleating. The 
question then is: does this multiplicity increase the probability of nucleation into a different phase? In other words, 
is there an 'entropic drive' for nucleation away from a state, leading the system to another state which, though 
individually not necessarily more favorable, belongs to an exponential family of states with free-energies comparable 
to the original one? The argument against an entropic drive goes as follows: While a system in a given phase is 
nucleating another phase A a , it needs to follow a certain number of steps. Why should this process be accelerated by 
the fact that there where many options of phases A±, A2 given that these themselves require different sets of steps 
to nucleate? To clarify the point, let consider a simple situation - a system at very low temperature, activating its 



escape out of the spherical crater V(r) of dimension d. in Fig. 13 Starting at the bottom, the particle follows (say) 



a Langevin process, and the probability evolves via a Fokker-Planck equation. 

P = V [TV + W] P (29) 




FIG. 13: Two trajectories escaping a crater. 
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The probability distribution was spherically symmetric at the start, and it will remain so, so we may go to spherical 
coordinates: 



P=' 



T ± ( r ^ dP \ +P d L U-x^C] + d ^ dP 
dr V dr J dr \ dr J dr dr 



Putting P = r d 'Pwe get the radial diffusion equation 

d 2 d 



P = 



T— + -(V(r)-T(d~l)lnr) 
dr z dr 



(30) 



(31) 



Thus, the system evolves in an effective potential V(r) — T{d — l)lnr. The second term in this effective potential 
derives from the multiplicity of excited states: (d — 1) lnr is the entropy of excited states, being the logarithm of the 
volume of a shell of radius r. We thus conclude that the multiplicity of excited states have the effect of lowering the 
barrier against nucleation. In particular, with very small T, but with Td = O(l) such that V(r) — T(d — l)lnr is 
monotonic decreasing, the particle will escape in finite time. 

An intuitive way of understanding this result is as follows: At low temperatures, for every trajectory that escapes, 
there are many unsuccessful trials which fall just short of escape, and return to the center. Although still extremely 
rare, trials reaching, say, 90% of the barrier height are exponentially more common than those which actually escape. 
Their frequency is, however, proportional to the surface of the 90% energy level. The situation for one out of an 
exponential number of metastable states is very similar. The nucleation of a droplet of of any other state of volume 
V will cost a free energy at most proportional to the surface. However, there are e vs ^ options of different droplets 
of free energy density /, and they will be accessed with total frequency J df e v ^~^^. These need not be grown 
simultaneously, but the many unsuccessful attempts in which a droplet grows and collapses will be this number of 
times as frequent as if there were only one other option. The existence of many paths (in the particle language) or 
metastable states (in the system language) does not make it easier to nucleate into any particular state, but it makes 
it easier to nucleate out of the initial deep well. In other words, if we lump the ensemble of other states as as a single 
entity (the liquid state, in fact), the probability of nucleating one of this ensemble (that is, of escaping a specific 
state) is just as if there were an extra term in the free energy of the liquid equal to —Tin \J df e VTl ^~ '^1 . Therefore, 
the usual nucleation argument tells us that this volume term will sooner or later overcome the surface term, and the 
system will leave the original state: this competition leads to escape only if there is an extensive complexity. 

Going back to the example of escaping a crater in large dimensions (such that T d is of order one), we may describe 
this at low temperatures either by an ensemble of extremal trajectories with different angles (the low temperature 



treatment of (29), or by a single escape trajectory for the radius (i.e. the low temperature limit as applied to Eq. 
(3D). The ensemble of escape routes at different angles in the original space, and the single escape trajectory in 
radial coordinates play respectively the role, played in an extended system by the nucleation of many phases and the 
nucleation of the 'liquid phase' that encompasses all of them. 

Technically, activation out of a state may be calculated by considering the extremal trajectory, starting in the state 
and moving away from it by nucleating another phase. The logarithm of the escape time is given by the action of the 
trajectory. In ordinary nucleation, the trajectory goes through a free energy maximum at the critical droplet radius, 
and the total action of the escape trajectory is of the order of the free energy barrier at the critical droplet. The 
situation here is different, because we have an exponential number of different trajectories moving out of the state, 
none of which needs to have a finite critical droplet radius. When we sum over all trajectories, it is the multiplicity of 
trajectories that decreases the action, in a manner entirely analogous to the way that multiplicity of states decreases 



the total free energy (cf Eq (28)). That is to say, we have a trajectory complexity that decreases the values of the 
dynamic action, and hence the logarithm of the escape time. If the action of an individual nucleation keeps growing 
with the surface of the droplet size, but the number of trajectories is proportional to the exponential of the droplet 
volume, then escaping will sooner or later win - even without a critical radius for an individual passage. 

We thus conclude that the state will be destabilized by the presence of exponentially many others, even if the 
passage to any individual one is extremely improbable. Note that this argument also implies that it is impossible to 
have an exponential number of states with super- Arrhenius lifetimes all the way down to zero temperature, because 
the Arrhenius escape probability will sooner or later dominate. 

The argument we outlined above is close in spirit to the mathematical results in [1] and to the construction in 
Ref. [4"0TPI2"] related to the point-to-set correlation function. What has been proven there, at least for lattice systems 
with discrete variables, is the following: if one considers a closed volume, and asks if a fixed configuration on the 
surface is correlated with the one in the center (a point-to-set correlation), the maximal radius r max at which this 
happens diverges if there is a divergent timescale. Translated into density profiles, this might seem to suggest that, for 
r max ~~ ^ o° the surface configuration determines the time-averaged configuration inside, and this would imply that the 
complexity is necessarily subextensive, since the number of different configurations of the boundary is proportional to 
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the exponential of an area. There is a problem with this argument, as the simple example of classical mechanics shows. 
In that case, giving the configuration at two times determines fully a trajectory in between, and yet, the number of 
possible trajectories may be exponential ( the Kolmogorov-Sinai entropy need not be zero). The reason for this is 
that, because of the Lyapunov instability, to know the configuration in the whole interval one needs to know those of 
the boundaries to greater and greater precision - the number of digits necessary being proportional to the interval of 
time considered times the Lyapunov exponent. In conclusion: the entropy of the bulk is not necessarily proportional 
to the boundary surface, even in the case that the bulk configuration is fully determined by the boundary one, but 
may be - in the case of exponential sensitivity to boundary conditions - also proportional to the enclosed volume. 



V. THE KAUZMANN/ RANDOM FIRST ORDER SCENARIO 

Assume that there is a temperature at which there is an equilibrium situation with permanent density modulations. 
We now ask: Are there other excited states with infinite lifetimes, that we may define in a time-averaged way? 
For example, a ferromagnet with periodic boundary conditions may be composed of two domains of positive and 
negative magnetization, separated by an essentially flat (at sufficiently high dimensions) interface. Another example 
is a crystal with a single straight dislocation. A third example was already mentioned above: in quasicrystals one may 
have isolated defects, the most famous being the 'decapod' defect in a Penrose tiling. In addition, in any non-periodic 
system, new states may be obtained from translations, if this does not violate the boundary conditions. We call these 
states because we are considering the ensemble of configurations visited in infinite times (the thermodynamic limit 
has been taken first) by a system that started with any one configuration with a given permanent defect. In other 
words, we are studying the problem modulo short-lived (such as particle-vacancy) excitations. 

Next, we ask about the free-energy of these states. In many cases, such as the macroscopic dislocations or domain 
walls mentioned above, their free-energy is infinitely (though non-extensively in the system's volume V) higher than the 
ground state. Their contribution to the equilibrium distribution hence vanishes. In other cases, it could conceivable 
happen, and the Wang-tile example cited above is an example, that an excited state of this kind has a free-energy 
difference of order one with respect to the ground state. In that case, the conclusion is that it will have a finite 
probability in equilibrium, and we are forced to conclude that the equilibrium distribution contains several 'ergodic 
components', mutually inaccessible (for a macroscopic system) in finite times. Note that the state with a single 
'serious' defect will typically have many short-lived activated 'benign' defects. 

We now ask about the number of these states, in terms of their free energy, at a given temperature. If we count how 
many of these states there are, regardless of their probability of appearing in an equilibrium sample, the logarithm of 
this number is the complexity E(T) (only weakly dependent on T). The probability of each state is proportional to 
e -^V(/ ,+/ < >) j £ Qr a sam pi e f volume V . The definition of E is extremely tricky: as soon as E(/) is of order one, the 
simple 'entropic drive' argument above says that the states cannot be stable. This is most clearly seen in the example 
of the Wang-tile system: as soon as we have a finite density of defects, and hence an extensive E, they are separated 
by finite distances and may be annihilated with local rearrangements. Strictly speaking, then, E is only well defined 
when it is subextensive, and it othe rwise counts states whose definition is timescale- dependent. 



As mentioned above (cf. eq. (28)), the whole question is whether 



f=Q T d f 



> 4. Only in the former 

f=o 



df 

case will a distribution frozen in the lower states be dominant. In the latter case entropy wins, states with higher 
free energy density dominate, and since these are exponential in number, they are necessarily unstable: we are in the 



liquid phase. Suppose then that 



< A holds, and we have a solid with permanent density modulations. As 

the temperature is increased, the system eventually becomes a liquid. How is the transition made? There are several 
possibilities: 



• Second order: the density modulations diminish with increasing temperature, and eventually vanish completely 
at a temperature T*. This means that the level qEA of the plateau in the autocorrelation decreases continuously 
and vanishes at a given temperature T*. This is what would happen in a ferromagnet, with the plateau value 
being the magnetization squared. 

• First order: a set of exponentially many (i.e. necessarily unstable) states suddenly begins to dominate the Gibbs 
distribution above a temperature T m . Such a transition is like melting. Note that if E suddenly becomes large, 
the states will necessarily be individually short-lived, and the system will jump to a liquid without high viscosity. 

• The E vs. / curve reaches zero with a finite derivative (3r (Figure [14] ). Because, as mentioned above, E is 
in general weakly dependent upon temperature, we must state that this has to happen at temperature Tr. At 



T > Tr the saddle point |? 



= rp begins to dominate. The saddle point complexity starts from zero and 

/=0 
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grows gradually, unlike the case of melting. Just above Tk, the states that dominate have small complexity 
and are, because of this, just barely unstable, the 'entropic pressure' being very small. This is the Kauzmann 
scenario, as revisited from the Random First Order perspective. 



A particular case of the above is when ^ 



timescale dependence, but only amorphous order at T = 0. 



= oo as T — >• 0. In that case, the system has super- Arrhenius 




FIG. 14: Complexity of states versus free energy density Y = f — /. Three situations within the Kauzmann-RFOT scenario. 
At small temperatures, the system is frozen in the lowest states (left panel). As soon as T > Tk an exponential number of 
metastable states dominate, their lifetime is necessarily large but finite, and the system becomes a liquid (center panel). Out 
of equilibrium, the system finds itself exploring states of free energy higher than those dominating equilibrium. If it does so 
in essentially a flat manner, we may assimilate the situation to a two-temperature scenario: T inside a state, and T e fj > T 
between states (right panel) 





constant (liquid) 


amorphous 




q<q c 




~~~~\ 





q>q c 



FIG. 15: The Renyi entropies of an ideal glass in equilibrium would be dominated by the lowest states (q > q c ), or the higher 
liquid- like states q < q c , with q c — x the Parisi parameter. For the supercooled liquid the situation is the same, except that 
q c > 1, so that typical states are liquid-like. 



A. A crucial hypothesis: from global to patch quantities. 

A standard result in statistical mechanics is that the entropy of a system is calculable, to leading order in the 
volume, as the limit of the patch (or 'block') entropy . This is because, even though the boundary conditions of each 
block are not those of the total system, one assumes that for short-range correlated systems the effect is of the order 
of the boundary surface, and hence subdominant in the large volume limit. If we assume that the same holds for 
the complexity, and furthermore, for all the Renyi complexities of finite q, then we may translate the predictions of 
theoretical scenarios for the entire probability distributions to observable statements concerning the distribution for 
patches. In other words, we are led to the assumption that the large-deviation function G — y\nV, expressed in 
terms of the normalized logarithmic frequencies , is the same, to leading order in V, as the one obtained for the 
whole system with volume V, and given boundary conditions. Note that this does not mean that individual profiles 
dominating patches of given volume V are the same as the state that would dominate a single patch with, say, periodic 
boundary conditions. 
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The Random First Order scenario is strictly speaking, at least for the moment, a mean- held notion, and any 
extrapolation to hnite dimensions is made at the price of phenomenological thinking (see, however, Ref.|27j). The 
argument above allows us to make the reasonable assumption that the predictions of RFOT for the number of states 
of the entire system may be applied to patches embedded in an infinite system as well. The RFOT theory predicts, 
at least for mean-field models, a series of metastable states with internal free-energy f a . For a given /, there is a 
quantity Af(f) of these states, with lnTV ~ VS sys (/), which defines the complexity S sys - a function that depends 
only weakly on temperature, because states retain their integrity when temperature changes. The assumption that 
all Renyi complexities of the whole system may be computed as a limit of the Renyi complexities for large patches 
means that we may identify T, sys and E. Thus: 



equil) 



e. 



-VP f a 




-Vqf a 



— qhi 



-VP fa 



V(q-1 

^ {in I df e n*U)-W] -qlnjdfe v ^U)-m 



V(q 



V(q-l) 



In / df e 



V[S(f)-Pq(f-f)] 



where 



/3/^--ln / ,//■, 1 ' ' ' 



Comparison with equation (27) yields, and we recognize q as the m parameter of replica theory 



(32) 



(33) 

Varying the Renyi 
For 
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parameter q is then exactly the procedure introduced by Monasson |29j to analyze these models (see Fig 
example, this means that if the system were a one-step replica symmetry breaking solution in the low temperature 
phase, with Parisi parameter x, then the Renyi complexities would be (see Fig 



14): 



K„ 







K q > 



for q < x 
for q > x 



(34) 



This is directly measurable if an equilibrium configuration is available (Fig [14]) . 

Experimentally, we thus have to proceed as follows. We first define an averaged profile, over a time r p of the order 
of r a . Note that r a is not the same for all patches, with the rarer ones expected to have shorter lifetimes. Hence, we 
must be careful of the fact that the number of rare states (i.e. the Renyi complexities for large q) may be dependent 
on the averaging time t p . 

Next, we either compute the Renyi complexities a la Grassberger-Procaccia, by scaling the results for different V 
and different e, or , alternatively, we set up a large deviation function plotting the logarithm of the probability of 
having a number of repetitions versus the logarithm of the number of repetitions, both normalized by V. 

A fingerprint of the RFOT scenario is then a non-zero crossover Renyi parameter q = x, implying that there are 
always rare patches that repeat exponentially infrequently. On approaching the transition from below x — > 1, and this 
means that non-repeating patches become more and more frequent, and eventually typical at T = Tk- Conversely, in 
the liquid phase just above Tr, most patches repeat exponentially rarely. However, there are a few ones that repeat 
frequently, but their number is so small that even their high repetition rate does not allow them to dominate the 
measure. 

Let us mention here that we may take the point of view of patches, rather than the whole system, for more 
complicated (and informative) large-deviation functions, depending on the patch overlap Q, defined, for example as: 



Qab = y J dr (p a (r) 



P)(Pb(r)-p) 



(35) 



One takes a patch at random, and counts the average complexity K\(Q) as in Equation (19), but restricting the sum 
to patches at distance Q. Finally, one averages over initial patches. Ki(Q) is the patch version of the Franz-Parisi 
effective potential [3D], a function closely related to the Parisi ansatz, for which there are also suggestions within the 
RFOT scenario. 
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VI. OUT OF EQUILIBRIUM: EFFECTIVE TEMPERATURES 



In practice, glasses are out of equilibrium, and the concept of temperature does not apply in principle to them. 
An old idea [31] has been to assume that the structure somehow "remembers" the situation of the liquid at the 
temperature it had before it fell out of equilibrium. The system would have one (or, in later versions, several) "Active 
temperatures" governing the slow evolution. On the other hand, when analytic solutions for the glassy dynamics of 
mean- field models became available [32] . it turned out that these systems possessed an effective temperature T e ft 
relating fluctuation and dissipation (of any observable) at low frequencies [33j . Indeed, T e ff would be measured by 
applying a low- frequency thermometer to the system out of equilibrium, an experiment that has been performed [34j . 
The interpretation for the appearance of this temperature is not yet fully clear, but it is not related to structures that 
are frozen in the form they had when the system crossed the transition, but rather to structures that form again and 
again. We know this because effective temperatures exist at times which are long compared to the autocorrelation 
time, so the system had more than enough time to restructure. 

A somewhat strong assumption is to suppose that this temperature arises because a system out of equilibrium 
explores metastable states with a probability at T e f / - a proposal that agrees with a hypothesis that had been made 
by Edwards for granular matter |35H37j - and that T e ff slowly evolves in a much longer (perhaps infinite) timescale 
towards the equilibrium value T. In this very strong form, the assumption is that if the system has metastable states 
with free energy f a , so that the equilibrium probability would be proportional to e~f"/ T ; out of equilibrium the 
probabilities would be those associated with a probability cx e~f a ' Te fi . In other words, if the probability of a certain 
configuration in equilibrium is pa qull \ the probability out of equilibrium should be: 



(out of equil) 
Pa 



(equil) 
ret. 



(36) 

Again, if we assume that what has been said for V can also be said for large V, then we should have: 

F(g-l)# g =ln^V«^M -glnf^V^-H (37) 

which means that, in every physical situation, the dependence on q may be written as 

iA ^K q = I3 eff [U(qf3 eff ) - U(f3 eff )} (38) 

where U(b) = i ^ Q e" b ^ a is a function that depends, upon this assumption, on the final parameters (temperature, 
pressure). All the history of the sample would manifest itself, under this hypothesis, only through the single parameter 



T e ff = l//3 e ff. In the limit q — > 1 (38) reduces to K\ — S|^ e// , the complexity evaluated in the point where £' = e ff- 
In order to check the hypothesis of an effective temperature of this kind, we should repeat different protocols, and 



for each determine the value of T e ff that will put the results in the form (38), with a single protocol-dependent T e // 
for all g's. This should be possible in principle, and if this is so, one could check that the patterns that dominate the 
measure at a given T e ff are the same, irrespective of the protocol by which it has been reached. 



VII. CONCLUSION 



In this paper we have discussed how to render observable several quantities from glass theory, in particular those 
that have their origin in mean-field like theories. We have discussed a correlation length that, we argue, should 
diverge for any solid, defined as a system in which, in spite of thermal fluctuations, there are permanent spatial 
density modulations. The construction has the strength of not being tied to any specific phenomenological scenario. 
Once we define states, order and lengths in a concrete and experimentally measurable manner, we are lead to the 
conclusion that the mere existence of a timescale that goes to infinity faster than an Arrhenius law implies, for finite, 
short-range interactions, that the lengthscale diverges. 

Defining further quantities associated with patches embedded within the system, we may answer the question of 
how, on the basis of the data from an ideal glass configuration, we could recognize if it corresponds to a Random First 
order/ Kauzmann scenario. Perhaps more relevant for realistic systems, we may also use the technique to decide, 
using the experimental out of equilibrium data, whether there is a flat exploration of metastable states, leading to an 
effective temperature. 
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